A new spinosaurid dinosaur species from the Early Cretaceous of Cinctorres (Spain)

A new spinosaurid genus and species is described based on the right maxilla and five caudal vertebrae of a single specimen from the Arcillas de Morella Formation (Early Cretaceous) at the locality of Cinctorres (Castellón, Spain). Protathlitis cinctorrensis gen. et sp. nov. is diagnosed by one autapomorphic feature as well as by a unique combination of characters. The autapomorphy includes a subcircular depression in the anterior corner of the antorbital fossa in the maxilla. The new Iberian species is recovered as a basal baryonychine. The recognition of Protathlitis cinctorrensis gen. et sp. nov. as the first baryonychine dinosaur species identified from the Arcillas de Morella Formation (late Barremian) from the same time as Vallibonavenatrix cani, the first spinosaurine dinosaur from the same formation in the Morella subbasin (Maestrat Basin, eastern Spain), indicates that the Iberian Peninsula was home to a highly diverse assemblage of medium-to-large bodied spinosaurid dinosaurs. It seems that spinosaurids appeared during the Early Cretaceous in Laurasia, with the two subfamilies occupying the western part of Europe during this period. Later, during the Barremian–Aptian, they migrated to Africa and Asia, where they would diversify. In Europe, baryonychines were dominant, while in Africa, spinosaurines were most abundant.

www.nature.com/scientificreports/ This formation, with a maximum thickness of 96 m, is composed of red clay, yellowish-white sandstone, grey marl, limestone and locally conglomerate. These lithofacies are arranged vertically and laterally, forming a set of depositional sequences limited by erosion surfaces. These sequences have a transgressive-regressive character and are composed of minor cycles of the parasequence type, which may present flood limits, subaerial exposure or erosion. The Morella Formation in the ANA fossil site has an average thickness of approximately 80 m, and a maximum of six depositional sequences have been identified in it. Following Bover-Arnal et al. 31 , the Arcillas de Morella Formation has an age of 0.745 M.a., although each of the six depositional sequences has an average age of 0.124 M.a. This is consistent with Milankovitch cycles of short orbital eccentricity of the Early Cretaceous astronomical periodicities 32 . Further analysis of the influence of astronomical parameters on sedimentology is in progress. All lithologies were deposited in flood plains, estuaries and beaches of a delta. The palaeocurrents come from the north and north-west from an emerged massif located at the centre of the Iberian Trough. This relief separates the Maestrat Basin from the Iberian-Levantine Basin. In the surroundings of Cinctorres where the fossil site is located, the thickness of the Arcillas de Morella Formation is 57 m.
The ANA site was discovered in 1998 by the geologist Ramón Ortí, but it remained unexcavated until 2002, when a palaeontological team formed by members of the Institut de Paleontologia Miquel Crusafont from Sabadell (Barcelona) and the Grup Guix from Vila-real (Castellón) unearthed the first fossil from the locality. The ANA site is located 37 m from the base of the Arcillas de Morella Formation and consists of approximately 2 m of grey-yellow sandy mudstones containing limonite and goethite crusts and nodules.
Systematic palaeontology. DINOSAURIA  Diagnosis: Baryonychine distinguished by the presence of the following unique autapomorphic feature: subcircular depression in the anterior corner of the antorbital fossa in the maxilla. In addition, Protathlitis cinctorrensis gen. et sp. nov. can be distinguished from other baryonychines on the basis of a unique combination of characters: in caudal vertebrae, the transverse process has only two fossae and one buttress, and they are posteriorly oriented; centra with a narrow ventral groove; circular outlines of the articular facets are clearly oval; prezygapophyses and postzygapophysis projected beyond the anterior and posterior rims of the centrum; and no type of element such as hyposphene below the postzygapophysis.
Type locality and type horizon: The ANA site is near Cinctorres (approximately 2 km), Castellón (Spain). The exposed beds belong to the Arcillas de Morella Formation, which is late Barremian (~ 127-126 Ma) 30,31 . Fossils were recovered in grey-yellow sandy mudstones containing limonite and goethite crusts and nodules that were deposited in a shallow low-energy estuarine environment (Figs. 1, 5). The universal transverse Mercator (UTM) coordinates of the ANA site are 30T YK 734714 4497540.
Description. Maxilla. 8ANA-109 is a nearly complete right maxilla (Figs. 2,3,4). It is broken into three separate pieces that fit together. It lacks the anterior part of the maxilla. The maxilla is anteroposteriorly elongate www.nature.com/scientificreports/ with a length of 430 mm preserved. The maximum width is 71 mm from the medial side of the medial shelf to the lateral wall and is 94 mm in height at its highest point. There are 16 alveoli separated by interdental plates. The maxilla lacks the anterior part, the zone of contact with the premaxilla. Part of the anteromedial process, the ascending ramus and the jugal ramus are preserved. Fragments of four teeth can be distinguished in the most posterior alveolus (Fig. 3). Additionally, some dentine laminae can be observed in other alveoli that form the lower part of the root of some teeth. The medial shelf is broken. The ventral margin is straight and tapers anteriorly with respect to the distal part in ventral view. The dorsal margin has two subcircular depressions, one situated in the anterior corner of the antorbital fossa (Fig. 3).
The lateral surface of the maxilla is clearly dorsoventrally convex. The lateral wall is pierced by numerous neurovascular foramina arranged in two rows. The lower one is very deep and bears three foramina with an elliptical outline. From each foramen, a groove runs anteroventrally. These grooves are straight or slightly parabolic and become deeper near the foramina. They penetrate the lateral wall of the ramus posterodorsally (Fig. 4). Among these foramina, a dorsally positioned row of two neurovascular foramina running anteroposteriorly can be distinguished. These holes are smaller. Additionally, they are elliptical in outline, a groove runs anteroventrally, becomes deeper near the foramina and penetrates the lateral wall of the ramus posteriorly. Both rows of foramina connect internally into a wide, oval-shaped canal (Figs. 3, 4). This canal runs anteroposteriorly and laterally to the length of the ramus and is parallel to the row of teeth. It originates on a lateral groove from the anterior medial wall of the antorbital fenestra. This canal is the promaxillary sinus 35 .
The maxillary alveoli extend far dorsally, but they are not totally vertical. They are oriented slightly lateral to the external wall of the maxilla. These alveoli are narrower anteriorly, and they become subquadrangular posteriorly. The interdental plates with the continuous interdental bone (sensu Currie 36 ) have an hourglass shape in ventral view. Four of the most posterior alveoli retain part of the root of the teeth, decreasing in size. The cross-section of these roots is elliptical. The last of the roots shows the lingual margin not preserved. It has been resorbed to let the replacement tooth migrate labially (Fig. 3), as is observed in other baryonychines 29 .
Vertebrae. The identification of the position of the isolated vertebrae in the vertebral column is difficult. In the following description, it is assumed that the vertebral column included forty-one caudal vertebrae or more. This www.nature.com/scientificreports/ is the general condition in basal tetanurans 37 . The centra of the vertebrae are fused to the neural arch, suggesting that the material belongs to a subadult or an adult form. The centra of the five anterior vertebrae are well preserved. They are spool-shaped with the side faces lightly compressed. They are amphicoelous, with the proximal face slightly more concave than the distal face. All of the recovered vertebrae are anterior to the "transition point", defined by the last caudal vertebra with transverse processes 1 and the first vertebra with elongated prezygapophysis 38 (usually occurring between Ca17-Ca19). The  Table 1), the total height of the centrum is slightly smaller than its length, and the centrum is longer than it is high. The shape of the articular faces is oval with the major axis in the dorsoventral direction. 4ANA43 shows a protuberance in the central part of the two articular faces of the centrum.
In lateral view, the distal ventral edge of the centrum is located below the proximal ventral edge, while the distal dorsal edge is rounded forwards to the junction with the neural arch. In the dorsal region of the centrum, just below the lateral processes, is the neurocentral suture, with a marked depression in 4ANA43 (Fig. 5). The more proximal vertebra has a deeper depression dorso-distally in both side faces, just below the neurocentral suture. This depression corresponds to the posterior region of a dorsal protrusion that is observed within the neural canal.
The distal articulation facets of the chevrons are more developed anteriorly than distal facets. There is a groove that divides both of them and continues anteriorly. The anterior facets are not as marked or differentiated as the distal facets. The ventral groove is shallower and narrower anteriorly. In the proximal region, there is a rounded and less developed crest that is located in the midline of the groove, which extends distally to approximately half the centrum.
The transverse processes are not completely preserved, except in 4ANA76 (Fig. 5), which maintains the left process, although it is deformed. The base of the processes is located slightly behind the midpoint of the vertebral centrum and is quite robust in section but flat and wide distally. The branch of the transverse process is directed slightly posterior and dorsally. Ventrally, the processes present a support ridge that is broad and slightly pronounced, which extends to half of their length.
The prezygapophysis of the anterior vertebrae is eroded. In 4ANA76, the facet joint of the right prezygapophysis is directed dorso-medially. In the prezygapophysis, there are well-developed spinoprezygapophyseal laminae. These laminae bind medially near the midpoint of the centrum, where the base of the neural spine arises. Among these laminae, there is a deep and narrow fossa. From the prezygapophysis, there are laterally prezygadiapophyseal laminae. They are bound to the proximal region of the transverse processes. The binding area of the process is located ahead of the union zone between the spinoprezygapophyseal laminae, forming a slightly concave lateral surface. This area is medially deeper near the laminae. www.nature.com/scientificreports/ The postzygapophyses of both vertebrae are eroded, and they do not allow observation of the facet joint. The basis of these is located in front of the distal edge of the centrum, and the angle between them is less than that between the prezygapophyses. At the centre of the postzygapophysis, a deep and narrow depression that extends dorsally is observed.
The neural spine is not complete in any of the vertebrae recovered. The base of the spine arises proximally at the junction of the spinoprezygapophyseal laminae, and caudally, it does not exceed the distal edge of the centrum.
The centra of the middle caudal vertebrae (3ANA83, 4ANA69, 5ANA78) ( Fig. 6; Table 1) have the same shape as the proximal caudal ones. The length of these vertebrae is greater than their height, as in the case of the proximal caudal ones, they are also amphicoelous, and the articular faces are oval with the dorso-ventral axis greater than the medial-lateral axis.
In lateral view, the distal ventral edge is not as low as in the proximal caudal vertebrae, and the depression of the neurocentral suture is not as marked. In these vertebrae, a deep depression was not observed in the 4ANA43 vertebra.  www.nature.com/scientificreports/ The articulation facets with the chevrons are more pronounced in the distal region of the vertebrae, and as in the case of the anterior vertebrae, and the facets are less marked in the proximal face. Ventrally, there is a slightly marked groove, but there is not a medial crest.
None of these vertebrae maintain complete transverse processes. The basis of the processes is long anteroposteriorly. The proximal edge arises near the anterior end of the centrum. The lamina of the process is thin in the dorsoventral direction, but in the 5ANA78 specimen, a less developed support crest is observed. The processes are directed posteriorly, and they are slightly dorsal.
The prezygapophysis is badly eroded, and the articular facet is not perceived. As in the case of the proximal vertebrae, the spinoprezygapophyseal laminae are highly developed, and they are joined together near the base of the prezygapophysis, forming a narrow but deep fossa anteriorly. In the 5ANA78 vertebra, the right prezygadiapophyseal lamina forms a thin ridge on the anterior edge of the process, as in the case of the anterior vertebrae forming a large lateral concave surface to the spinoprezygapophyseal lamina.
The postzygapophyses are long and protrude slightly behind the distal edge of the centrum. The articulation facets are large, and they are lateroventrally oriented at an angle similar to that of the anterior vertebra. Dorsally to the postzygapophysis, a deep and elongated depression is observed that is directed towards the neural spine.
The neural spine is not complete in any of the medial vertebrae. The base is slightly shorter in the anterior-posterior direction than in the case of the anterior vertebrae, and it originates from the joining of the spinoprezygapophyseal lamina, which is anterior to the midline of the centrum.
Phylogenetic results. Using the phylogenetic analysis program TNT 39 , we performed a "traditional search" with tree bisection reconnection. Wagner trees with a random seed of 1 and 9,999 replicates with 10 trees saved per replication were used. All characters were equally weighted and treated as unordered. TNT revealed 9,386,515 rearrangements. As a result, one most parsimonious tree (MPT) was recovered (Fig. 7) with a length of 157 steps, a consistency index (CI) of 0.8089, a homoplasy index (HI) of 0.1911, and a retention index (RI) of 0.8598 (the 'describetree' command in PAUP 4.0a build 167 40 ).
Mateus and Estraviz-López 1 described a new Spinosauridae dinosaur, Iberospinus natarioi. They pointed out that the fossil shows some characteristics of baryonychines, suggesting a closer relation to the group but placing it outside of both the Baryonychinae and the Spinosaurinae. This suggests that Iberospinus would be a sister taxon to Camarillasaurus (outside both subfamilies), or perhaps both taxa are the same (synonymy). Future research should focus on this topic.

Discussion
Maxilla. 8ANA-109 shares with Spinosauridae the proportionally larger medial shelf (anteroposteriorly and mediolaterally) compared with the rest of tetanurans 29 . Additionally, it shares with this family the posterior decrease in the size of the alveoli, the conidont dentition and the veined/anastomosed enamel texture 29 .
Sereno et al. 14 recognized some features that differentiate Baryonychinae from Spinosaurinae. The former has curved tooth crowns, teeth with fine serrations, external naris retracted to the first half of the maxillary tooth row and an elevated number of maxillary teeth (22 in Suchomimus). In contrast, Spinosaurinae has unserrated teeth; tooth crowns slightly curved or straight; external naris retracted farther caudally and fewer maxillary teeth (12 well-spaced maxillary teeth in Spinosaurus (MSNM V4047) and 11 in Irritator).
The elevated number of maxillary teeth (almost 16), curved crown and teeth with fine serrations situate 8ANA-109 as a member of Baryonychinae. In this subfamily, there are four valid genera, Baryonyx, Ceratosuchops, Iberospinus and Suchomimus 16 . Only Baryonyx 3 and Suchomimus 16 have maxillary remains. The 8ANA-109 fossil differs from Suchomimus 16 because it has a subcircular depression in the anterior corner of the antorbital fossa that is absent in the African dinosaur. This character in not preserved in Baryonyx 3 .
The shape of the alveoli is subcircular, as in other spinosaurids, but differs from the more elliptical alveoli in a fossil maxilla (CPI 477) from Igea (Spain) 29 .
8ANA-109 shows two replacement teeth at different stages of development within the same alveolus, as is seen in other theropods, such as Irritator, Oxalaia, Iberospinus and the CPI 477 maxilla 29 .
Vertebrae. The vertebral centra have the typical spool shape that is present in most theropods, with a U-shaped cross section, unlike in the case of Baryonyx walkeri, where the caudal vertebral centra are highly compressed 3 . Additionally, in the ANA theropod specimen, there is no type of element such as hyposphene below the postzygapophysis, as in B. walkeri 3  www.nature.com/scientificreports/ Caudal vertebrae from the ANA site differ from those of Riparovenator milnerae 2 in the shape of the neural arch. In R. milnerae, transverse processes start from the prezygapophysis, widely occupy the length of the centrum, and maintain their width transversally. Additionally, the prezygapophysis is elevated from the centra, the postzygapophysis is more developed posteriorly, and at the base of the transverse processes, the centrodiapophyseal lamina is absent in the anterior neural arches.
Additionally, neural arches of Spinosaurus 41 are very different from those from the specimen from the ANA site, with a notably complex array of vertebral laminae and fossae in the proximal caudal vertebrae. In Spinosaurus, the transverse process has three centrodiapophyseal fossae and two centrodiapophyseal laminae, but in the ANA specimen, the transverse process has only two fossae and one buttress. Additionally, the anterior caudal vertebrae of Spinosaurus exhibit wide and shallow ventral grooves that differ from the centra of the ANA specimen, which have a narrow ventral groove. Additionally, the transverse processes in Spinosaurus are projected perpendicular to the centra, while in P. cinctorrensis, they are posteriorly oriented.
Ichthyovenator 13 caudal vertebrae are clearly different. They have robust and elongate transverse processes that extend posterolaterally and dorsally. The first transverse process is sigmoid in shape in dorsal view. Prezygodiapophyseal and anterior centrodiapophyseal laminae are well marked and define a deep prezygapophyseal centrodiapophyseal fossa unknown in other theropods 13 .
The caudal vertebrae of Siamosaurus 15 (Phuwiang spinosaurid B), Camarillasaurus 4 and Iberospinus 1 are all similar 15 . These specimens mainly differ from Protathlitis in that they have two buttresses and three fossae below the transverse process of the caudal vertebrae. V. cani 5 is different in size, and ANA fossils are at least twice as large as from Vallibona. Another difference is the circular outlines of the articular facets that are approximately as wide as they are high in V. cani 5 , while in a similar position (Ca9), those of P. cinctorrensis (3ANA83, 4ANA69, 5ANA78) are clearly oval. The same occurs with the neural canal. The anterior and posterior edges are at the same level in Vallibona specimens, while in ANA vertebrae, the posterior edge is lower than the anterior edge. Postzygapophyses are short, projecting only slightly beyond the level of the posterior articular facet, in contrast with those from ANA, which clearly surpass the posterior articular facet. Finally, the transverse process of P. cinctorrensis is longer.
Suchomimus tenerensis 14,16 vertebrae are similar to those of Protathlitis, but differ because the prezygapophysis is more dorsally projected and the postzygapophysis is shorter (slightly projecting beyond the posterior rim of the centrum) than in P. cinctorrensis.
The ventral surface in the anterior caudal vertebrae has a groove in ANA fossils, as in Camarillasaurus 4 , Vallibonavenatrix 5 and Spinosaurus 41 . This differs in Baryonyx 3 , Suchomimus 14 and Ichthyovenator 13 with a flat ventral surface.

European evolution
Spinosaurids have been recovered in the western part of Europe, mainly in Portugal, Spain and the United Kingdom 1,2,5,22,29 . However, the most common fossils are teeth 23 23 and Barker et al. 2 , among others. It also seems clear that in the western part of Europe during the Barremian, at least two different taxa of spinosaurids coexisted 23 . All the findings produced to date seem to indicate that the group originated in western Europe during the Late Jurassic or Early Cretaceous and then expanded to Africa and Asia. For example, in the Iberian Peninsula, spinosaurids are still the most common theropods in the Lower Cretaceous formations (from the upper Hauterivian to Aptian 27 ). Another example is that both baryonychines and spinosaurines are documented in different areas, such as the Cameros and Maestrat basins 23 during the Early Cretaceous (Barremian-Aptian). During the Barremian-Aptian, spinosaurids moved from Europe to Africa and diversified there during the Aptian-Cenomanian period. Several migratory routes have been suggested between Gondwana and Laurasia during the Early Cretaceous, such as the ' Apulian route' 44,45 or during the Late Cretaceous 46 based on ichnites and skeletal remains.
Therefore, a possible scenario would be the following: spinosaurids appeared during the Early Cretaceous in Europe and diversified, appearing as members of the two subfamilies and occupying a large part of western Europe. Later, during the Barremian-Aptian, they migrated to Africa and Asia, where they diversified again. In Europe, baryonychines were dominant, while in Africa, spinosaurines were the most abundant.
Finally, previous works suggested that spinosaurids inhabited inland areas far from the coast or freshwater environments with marine influence or next to coastal zones 23,27,47 . In the case of ANA dinosaurs, they inhabited areas next to coastal zones.
These fossils belong to Spinosauridae because they have a larger medial shelf, alveoli decreasing in size posteriorly, conidont dentition and veined/anastomosed enamel texture. Additionally, they belong to Baryonychinae because they have curved tooth crowns, teeth with fine serrations, external naris retracted to the first half of the maxillary tooth row and an elevated number of maxillary teeth. Inside this subfamily, one autapomorphy supports the validity of the type species (subcircular depression in the anterior corner of the antorbital fossa in the maxilla). In addition, Protathlitis cinctorrensis gen. et sp. nov. can be distinguished from other baryonychines on the basis of a unique combination of characters: in caudal vertebrae, the transverse process has only two fossae and one buttress, and they are posteriorly oriented; centra with a narrow ventral groove; circular outlines of the articular facets are clearly oval; prezygapophyses and postzygapophysis projected beyond the anterior and posterior rims of the centrum; and no type of element such as hyposphene below the postzygapophysis. Regardless of the phylogenetic analysis carried out, P. cinctorrensis is clearly nested as a basal baryonychine (Fig. 7).
The establishment of this new European species seems to confirm that spinosaurids appeared during the Early Cretaceous in Laurasia, with the two subfamilies occupying the western part of Europe. Later, during the www.nature.com/scientificreports/ Barremian-Aptian, they migrated to Africa and Asia, where they would diversify. In Europe, baryonychines were dominant, while in Africa, spinosaurines were the most abundant.

Methods
Phylogenetic analysis. A phylogenetic analysis was conducted using a modified version of the Sereno et al. 16 dataset (Supplementary File S1) based on a 120-character list (49% cranial, 51% postcranial). We revised the dataset in its matrix, and we added one missing vertebral character to the specimen Baryonyx walkeri 3 (character 82: Caudal vertebrae, anterior, morphology of ventral surface changed from (?) to flat (0)). The matrix includes 15 × 148 operational taxonomic units (OTUs). For Protathlitis cinctorrensis gen. et sp. nov., twenty-five (15 cranial and 10 postcranial) characters were scored from the available information from the specimen (holotype). This modified matrix was analysed using TNT v.1.5 software 39 . All characters were equally weighted and treated as unordered. Using TNT software, we calculated the 'Bremer Support' , and we identified which synapomorphies supported each node. Only unambiguous synapomorphies were calculated, and branches with no possible support were collapsed ("rule 3"). This makes the results more conservative. Ceratosaurus, as a tetanuran outgroup, and Allosaurus, as a non-spinosauroid theropod, were employed as constraints.

Data availability
All data generated or analysed during this study are included in this published article and its supplementary information files.